Decay of polarons and molecules in a strongly polarized Fermi gas 
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The ground state of an impurity immersed in a Fermi sea changes from a polaron to a molecule 
as the interaction strength is increased. We show here that the coupling between these two states 
is strongly suppressed due to a combination of phase space effects and Fermi statistics, and that it 
vanishes much faster than the energy difference between the two states, thereby confirming the first 
order nature of the polaron-molecule transition. In the regime where each state is metastable, we 
find quasiparticle lifetimes which are much longer than what is expected for a usual Fermi liquid. 
Our analysis indicates that the decay rates are sufficiently slow to be experimentally observable. 



PACS numbers: 03.75.Ss, 05.30.Fk, 67.85. Lm 

The concept of a quasiparticle, developed by Landau 
in his theory of Fermi liquids [1], is fundamental to our 
understanding of low-energy excitations in many-body 
systems It has been successfully applied to describe 
systems spanning many energy scales ranging from elec- 
trons in a metal to atomic nuclei and quark-gluon plas- 
mas. In a recent experiment on ultracold gases, quasipar- 
ticle physics has been analyzed in depth by studying the 
properties of a few impurities immersed in a Fermi sea 
(FS) Theoretically, one has presently a good under- 
standing of the weak coupling limit of this problem where 
the ground state is a screened impurity, a state which is 
referred to as a polaron, and of the strong coupling limit 
where the ground state consists of a molecule formed by 
the impurity and a particle in the FS [H [S] . In the inter- 
mediate regime, the ground state changes from a polaron 
to a molecule as the interaction is increased, but the na- 
ture of this transition which involves a reorganization of 
the FS is not clear at present. The polaron-molecule 
transition is closely related to the problem of mixtures of 
dilute protons in neutrons at subnuclear density, which 
are encountered for example in stellar collapse. The pro- 
tons may form either polarons or two-body bound states 
(deuterons) with the majority neutrons, in analogy with 
what is considered in the present paper. 

We address here the polaron-molecule transition us- 
ing a diagrammatic expansion in the number of holes 
in the FS. The coupling between the polaron and the 
molecule is shown to scale with their energy difference 
Aw as |Aw|^/^. The high power is a result of phase space 
effects and Fermi statistics. It should be compared with 
the usual Ao;^ scaling of the quasiparticle damping rate 
in a Fermi liquid This scaling implies long lifetimes 
and a first order transition between the polaron and the 
molecule. Our analytical insights are confirmed by nu- 
merical calculations, and in the conclusions we discuss 
how the decay can be measured. 




FIG. 1. Polaron to molecule (a) and molecule to polaron (b) 
decay processes in presence of a fully-polarized Fermi sea. 



We consider a single J, impurity (which may be bosonic 
or fermionic) with mass immersed in a FS of f 
fermions with mass and density = kp/Gir'^. The 
interaction between the and f particles is short-ranged 
and described by the scattering length a, whereas the 
interactions between identical fermions can be ignored. 
Various approaches locate the polaron-molecule transi- 
tion at the critical coupling l/kpUc ~ 0.9 for — 
mi m [SHS]. In the strong coupling regime 1/kpa > 
l/kpUc, a polaron with energy ujp (we take /i = 1) is un- 
stable and will decay into a molecule with energy ujm by 
removing a particle from the FS. Due to the conservation 
of both momentum and energy, the decay must involve 
the creation of an additional particle-hole pair [i]. The 
leading decay channel is therefore the three-body process 
illustrated in Fig. [I] (a). Likewise for 1/kpa < 1/kpac, 
the molecule decays into a polaron by adding a fermion 
to the FS and creating an additional particle-hole exci- 
tation. This is again a three-body process as illustrated 
in Fig. 0(b). 

A polaron with momentum p is described by the imag- 
inary time propagator G;(p,t) = -(TV[ap;(r)a^^(0)]). 

Here, creates a J, particle with momentum p, Tr de- 
notes time ordering with respect to r G [0; 1/T[ with T 
the temperature (we take kg = 1) and (. . .) the thermal 
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FIG. 2. (a) The one-hole polaron self energy Ep-*. (b) 

(2) 

The two-hole polaron self energy Sp . Thin solid lines in- 
dicate G°, dashed lines the interaction, double wavy lines the 
molecule D, • the molecule-atom coupling strength g, and 
squares the scattering matrix T2. 



average. Its Fourier transform reads 

G;(p,2)-i = G°(p,2)-i-I]p(p,z) 



(1) 



where G'5J(p, z) = 1/iz - ^pi), = /Inia - iia, and 
/icr is the chemical potential for a =t,4. species. We take 
the ideal gas value — ep = kp/2m^ whereas /z^ plays 
no role as there is no macroscopic population of the J, 
states. The self energy Ep describes the effects of inter- 
actions. To obtain a systematic description of the decay, 
we expand the self energy as 



I]p(p,z) = EW(p,z) + S^^^(p,z) + ... . 

Here, Ep ■* involves processes with n holes in the FS. 

The term Sp'' shown in Fig. [2] (a) corresponds to the 
commonly used ladder approximation. The energy of a 
zero-momentum polaron obtained from 



(2), 



(2) 



wp = Sp''(p ~ 0,ujp + iO-i 



(3) 



(we drop the infinitesimal iO+ in the following) is essen- 
tially indistinguishable from the Monte-Carlo result [31|3]- 
[T^ . However, ([s]) yields a real energy for any cou- 
pling strength. To obtain the damping in the metastable 
regime it is therefore necessary to include processes in- 
volving creation of at least two holes. 

The Feynman diagrams for Ep' corresponding to the 
decay process in Fig. [T] (a) are given in Fig. [2] (b) . They 
describe the scattering of a p = polaron creating two 
holes with momenta q and q', and leaving one fermion 
above the FS with momentum k. The remaining fj, pair 
can then form a molecule with momentum q -I- q' — k. 

The propagator of the open-channel, Feshbach 
molecule is D{q,T) = -(r^[6p(T)6|,(0)]) with bl = 

/rf^#<7ap/2+q;ap/2-qt ^nd d^q = d^q/{2T:f. Ignoring 
finite range effects, its wave function in vacuum reads 
+ q^a^). The Fourier transform of the 
molecule propagator can be written as 

D{p, z)'^ = i?o(p, - Em(p, z) 
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FIG. 3. (a) The molecule propagator in the ladder approxi- 
mation, (b) Diagrams corresponding to the leading molecule 
decay channel. Wavy lines denote the molecule in the lad- 
der approximation, dotted lines with the operator 6 or b^, 
and double lines the polaron propagator (G4. — GJ). 



where Tim is the molecule self energy and 

,2 T^iv.z) 



(5) 



describes the propagation of a fj, pair in the ladder 
approximation. It is illustrated in Fig. [s] (a). Here 
f{x) — [exp(a;/r) + 1]^^ is the Fermi function and 
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l-/(ep-qT) 



(6) 



is the atom-molecule coupling. The object entering the 
diagrammatic analysis is ^(p, 2;)^-D(p, 2), which describes 
the scattering of atoms mediated by the Feshbach 
molecule |13) . In the following we use the vacuum limit 
of ([g]), g = —\/2'K /m^a with to,- = to-^to^/(to^ -|- to-|-). 

Using the fact that there is no macroscopic population 
of the 4, state and that the f fermions form an ideal gas, 
the frequency sums in the two diagrams can be performed 
analytically by contour integration. This amounts to 
evaluating frequencies at the on-shell f energies. The po- 
laron decay rate is given by Fp = — ImEp'^ (0, ujp). When 
< Aw = ojp-ojM < ^F, the Gj's andTa's inFig.[2](b) 
are off resonant and therefore real. The only contribu- 
tion to the imaginary part giving rise to damping comes 
from the Feshbach molecule, which in the vicinity of its 
pole can be described by 



Z 



M 



(7) 



Here, Zm and m\j are the molecule residue and effective 
mass respectively. The contribution to Fp from the two 
diagrams in Fig. [2] (b) at T = can after some algebra 
be written in the symmetric form 

rp = ^3^ ( d^qd^kd^q'[F{(i,k,u;p)- F{q',k,u;p)f 



xs {Aoj + + ^q't - ^fct - (q + q' - k)V2? 



(8) 



We have defined F(q, k, w) = T2 (q, uj + Cgt)G? (q - k. 



(4) ~ Cfct)- Here and in the following, we take q,q' < kp 
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and k > kp. The summation over two holes and one 
particle in (|8| can be interpreted as a reorganization of 
the FS to contain one less particle. 

For Aw <^ ep, the holes and particles involved in 
the scattering process illustrated in Fig. [T]^a) are located 
around the Fermi surface, \.e. q c:^ k ~ k' ~ kp. Also, 
the molecule momentum is small since /2m\j < Aw. 
Using this and ignoring for the time being the matrix 
element in (|8|, we get the integral 



P 



2m 



M 



25/2 

105^ 



(9) 



Thus, the available phase space of the three-body process 
gives a factor Auj'^^^ for the decay rate: the (5-function 
removes one energy integral in ([9|, a factor Auj comes 
from each of the two remaining degenerate energies, and 
a factor Aoj'^/^ comes from the molecule. Let us now 
focus on the matrix element F(q, k,wp) — F(q',k,wp) 
in ([s]). The momenta q, q', and k form an equilateral 
triangle for Aw — and the matrix element vanishes, 
since F(q, k, w) only depends on the angle between q 
and k. For Aw <^ ep, we can expand the matrix element 
in the deviations of the triangle formed by q, q', and k 
away from the equilateral shape. Using this in (|8|) yields 



Tp ^ Zjsfkpa 



Aw 



9/2 



ep. 



(10) 



Thus, the Fermi antisymmetry when swapping the hole 
momenta q and q' gives an additional factor Aw to the 
decay rate. The resulting polaron lifetime l/2rp diverges 
much faster than 1 / Aw close to the transition point kpUc- 
From ( 10 1 it follows that the matrix element between 



equal energy states which contain a molecule or a polaron 
vanishes at the transition point. This means that there is 
no avoided crossing and that the transition between the 
two states is of first order, as was suggested previously g]. 

Consider now the regime Aw < {1/kpa < 1/kpac)- 
Here the molecule decays into a polaron via the three- 
body process illustrated in Fig.[l|b). The corresponding 
Feynman diagrams are shown inTSTb) . After the molecule 
has split, the impurity scatters with the FS leaving one 
hole with momentum q and two particles with momenta 
k and k'. The impurity can then form a polaron with 
momentum p = q — k — k'. The decay rate for a zero 
momentum molecule is given by Tm = — ImSM(0, wm), 
and calculations like the ones for the polaron show that 
for r = and |Aw| ^ cp the two diagrams give 

Tm = ^ f Skd'k'd^ [C(q,k,WM) - C(q,k',WM)]' 



C(5 (|Aw| -f - a-T - 6't - (q - k - k')V2m^) 



(11) 
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FIG. 4. Polaron— ^molecule and molecule— >polaron decay 
rates Vp and Tm as a function of the energy difference |Aw| 
for m-t- = mi and T = 0. The thin dashed hues are the 
predicted asymptotic behaviours, Vp oc (mX/)'^'^^ALi;^''^ and 
Tm oc (mp)^''^ Aw^''^. Inset: Aw as a function of l/kpa. 



Here, C(q, k, w) = G^k, w - ^kt^ik - q, w - ^fet + ^^t) 
and Zp is the polaron residue. Like for the polaron decay, 
one can derive from ITll) that 



M 



Zpkpa 



I Aw 



9/2 



(12) 



for I Aw I ^ ep. Again, this means that the molecule is 
long lived close to the cross-over. 

In Fig. Q we show calculations of the decay rates ob- 
tained from Eqs. (|8|) and (11) for the — case. In 



the numerics we have used 



T2(P,W) 



2na/mr 



1 — y/2mra'^{p'^ /2mi^i — w — ep + gsn^) 



which describes molecules with energy 



ujm = ^fi^/ (2mra^) — ep -t- gsn^. 



(13) 



(14) 



Here 



.93 



27r/i^ 03/7713, ma = m^{m^ + m^)(2m^ -I- 
m^), and 03 = 1.18a is the molecule-atom scattering 
length [TU [TS] . Equation ( 14 ) is exact in the limit 



1/kpa — ^ 00, and accurate all the way to the critical 
coupling 1/kpac for the case = [4j. The po- 
laron energy is obtained from ^ , while m*p and are 
taken from Ref. [3]. The numerical results shown in Fig.|4] 
match nicely with (10) and (fl2|) for |Aw| <^ ep. In the 



fits, we have taken into account the factors (mp)"^/^ and 
(m^)^/^ coming from the (5-functions [see For larger 
I Aw I, there is as expected some discrepancy between the 



numerical results and (10) and (12) since the momenta 



involved in the scattering processes are no longer close to 
the Fermi surface. 
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Our analysis shows that the decay of polarons and 
molecules will happen on timescales of order 10 — 100ms. 
To observe the decay, one could for instance sweep 
the scattering length across Oc- Zener type of argu- 
ments predict that the polarons/molecules will survive as 
metastable states after the sweep (TH]. Radio-frequency 
spectroscopy [3] or optical probes can then be used 
to map the exponential decrease in time of the po- 
laron/molecule population. Alternatively, one can ex- 
tract the effective mass from collective mode frequen- 
cies [T7] . This technique has been successfully applied at 
unitarity [TH| . By measuring the change in the frequency 
as a function of time, one may extract the decay rate of 
the metastable state. 

The polaron can decay to deeply bound closed chan- 
nel molecules via three-body processes analogous to the 
ones considered here for the decay into the weakly bound 
Feshbach molecule. This decay is however strongly sup- 
pressed because the deeply bound molecules are spatially 
much smaller than the Feshbach molecule, whose size is 
^ a. Since the decay requires two t fermions to be at 
distance of order of the molecule size, the Pauli principle 
suppresses decay into deep molecular states much more 
efficiently than into the larger Feshbach state [TH]. The 
Feshbach molecule furthermore has a large component 
in the open channel for broad resonances, which justifies 
our use of a single channel theory in the present paper. 
Contrary to this, the deeply bound molecules are mostly 
in the closed channel where the fermions are in different 
spin states. This should suppress the decay even further 
as it must involve spin flips. Indeed, no decay to deeply 
bound states has been observed for a strongly polarized 
Fermi gas [3] , and very long lifetimes of order tens of sec- 
onds have been measured for balanced two component 
Fermi gases [50]. Contrary to this, our results show that 
the decay into the weakly bound Feshbach molecule is in 
the 10- 100ms range, confirming that it is the dominant 
loss process. 

One can include more interaction effects in the decay 
for instance by replacing GjJ with in the matrix el- 
ements or by introducing additional scattering events. 
This will change the quantitative value of the matrix el- 
ements for the scattering processes but will not change 
the F ~ |Aa;|^/^ scaling as this comes from the combined 
effects of kinematics and Fermi statistics; it can in fact 
be obtained from a Golden Rule calculation. 

For = mi and T = 0, Monte-Carlo calculations 
predict at small impurity densities a phase separated 
state before polarons may decay into molecules 12J. We 
expect however that non-zero temperatures will stabilize 
the polaron state against phase separation due to the en- 
tropy of mixing. Moreover, the boundary between mixed 
and separated states depends on the ratio m^/m^ |21j . 
By appropriately selecting the atomic species, one may 
bring the critical crossing point back to a physically ob- 
servable region. Alternatively, one could use bosonic im- 



purities, for which stability against phase separation has 
recently been predicted [22]. It would therefore be in- 
teresting to extend experiments with large population 
imbalances to mass-imbalanced fermionic mixtures, or to 
Bose-Fermi mixtures. 

In conclusion, we have considered the leading three- 
body processes involved in the polaron-molecule cou- 
pling. Our analysis shows that the coupling is strongly 
suppressed due to a combination of phase-space and 
Pauli blocking effects, vanishing as a power law close 
to the transition. This yields very long lifetimes of the 
metastable states, and implies a first order transition be- 
tween the two ground states. Our results suggest new 
directions for experiments with polarized atomic gases. 
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